library(tidyverse)
library(sf)
library(terra)
library(openxlsx)
library(mapview)
library(RColorBrewer)
library(osrm)
library(movecost)
library(fs)
library(stringr)
library(synthdid)
mapview::mapviewOptions(fgb=FALSE)
load("data/china.RData")
日本の1992年と2018年における夜間光
ntl1992<-terra::rast(x="data/9828827/Harmonized_DN_NTL_1992_calDMSP.tif")
ntl2018<-terra::rast(x="data/9828827/Harmonized_DN_NTL_2018_simVIIRS.tif")
ntl1992 <- ntl1992 %>%
terra::crop(y=ext(c(123.1,146.9,24.0,45.7)))
ntl2018 <- ntl2018 %>%
terra::crop(y=ext(c(123.1,146.9,24.0,45.7)))1992年
terra::plot(ntl1992)
2018年
terra::plot(ntl2018)
\[ \sum \]
jpn1<-sf::read_sf(dsn="data/JPN_shp/gadm41_JPN_2.shp")
mapview::mapview(x=jpn1,
zcol="NAME_2",
color="transparent")